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We obtained the spectral function of the graphite H point using high resolution angle resolved 
photoelectron spectroscopy (ARPES). The extracted width of the spectral function (inverse of the 
photo-hole lifetime) near the H point is approximately proportional to the energy as expected from 
the linearly increasing density of states (DOS) near the Fermi energy. This is well accounted by 
our electron-phonon coupling theory considering the peculiar electronic DOS near the Fermi level. 
And we also investigated the temperature dependence of the peak widths both experimentally and 
theoretically. The upper bound for the electron-phonon coupling parameter is 0.23, nearly the same 
value as previously reported at the K point. Our analysis of temperature dependent ARPES data 
at K shows that the energy of phonon mode of graphite has much higher energy scale than 125K 
which is dominant in electron-phonon coupling. 

PACS numbers: 81.05.Uw, 63.20.kk, 73.20.At, 79.60.-i 



I. INTRODUCTION 

Fermi liquid theory^ (FLT) is thought to be one of 
the most successful theories for describing the behaviors 
of electrons in solids, especially electrons near the Fermi 
energy in metals at low temperature. The success of the 
FLT in metallic systems naturally raises an issue on how 
far the FLT scheme can be applied to other condensed 
matter systems. Related to this question, there is a long- 
standing controversy on whether electrons in graphite, 
a 2 dimensional semi-metallic system, can be described 
within the FLT scheme or not. According to FLT, the 
lifetime of an electron due to electron-electron interac- 
tions is inversely proportional to the square of the bind- 
ing energy. Therefore, measurement of the lifetime as a 
function of the binding energy of an electron would be a 
direct test of the validity of FTL in graphite. 

Experimental results do not seem to show evidence for 
Fermi liquid behavior of electrons in graphite^. In fact, 
the inverse lifetime measured by 2 photon photoemission 
experiments (2PPE) conducted on natural single crystal 
graphite and highly oriented pyrolytic graphite (HOPG) 
appears to increase linearly as a function of the binding 
energy^. The observed peculiar behavior in the energy 
dependence of the inverse lifetime was discussed in terms 
of the peculiar dispersion of plasmon^ or electron-electron 
interaction in combination with the band structure of 
graphite-5^. 



However, electron-phonon coupling (EPC), one of the 
most fundamental interactions in solids, has not been 
considered in the discussion. On the theoretical side, 
very little work can be found on the EPC in semi- metals 
even though it has been well developed and widely stud- 
ied in metallic systems^. Only very recently has some 
theoretical models for graphene appeared^i^i^. On the ex- 
perimental side, electron lifetime was measured only for 
the energies larger than the maximum phonon energy of 
graphite (~ 200 meV)^^^ in the 2PPE experiments^!^. 
Therefore, to address the lifetime issue due to EPC in 
graphite, one may need two requirements. First, the ex- 
perimental data must show the lifetime of quasi-particles 
sufficiently close to the Fermi energy, less than the maxi- 
mum phonon energy of graphite. Second, a proper model 
that considers the electron-phonon interaction contribu- 
tion to the quasi-particle decay should be developed. In 
regards to the second point, models developed for metals 
have been used in the analysis of angle resolved photoe- 
mission (ARPES) data on graphite due to the lack of 
theoretical EPC models for semi-metals 

To address the issue of the quasi-particle dynamics and 
EPC in graphite, we performed high resolution ARPES 
experiments on high quality natural single crystal and de- 
veloped a theoretical model that considers linear density 
of states (DOS) near the Fermi energy^^. Our previous 
work was performed near the i^-point and showed a rel- 
atively small EPC constant A = 0.20. To extend our 



previous we have obtained high resolution ARPES data 
from the H point to determine if it also has a small EPC 
constant. In addition, we also performed temperature 
dependent studies near the K point. The temperature 
dependent data is compared with a theoretical model 
that fully considers the graphite DOS. Properly extracted 
peak widths are well understood within our EPC model 
with a linear DOS near the Fermi energy and shows a 
small EPC constant of less than 0.23. 



II. THEORY 

We first consider the theoretical side of the quasi- 
particle dynamics in graphite. In this section, we discuss 
possible decay channels for quasi-particles in graphite. 
First, it will be discussed how the lifetime of quasi- 
particles in graphite can be affected by EPC. We will 
formulate the self energy of quasi-particles based on the 
linear DOS of graphite, for zero temperature in section 
1 and for a finite temperature in section 2. The latter is 
to establish the foundation for estimating the EPC con- 
stant through temperature dependent studies. Second, 
we will discuss other decay channels such as electron- 
electron scattering, electron-plasmon scattering, and im- 
purity and defect scattering. Through these discussions, 
we wish to establish that the dominant scattering mecha- 
nism for quasi-particles in graphite comes from the EPC. 



A. Electron-phonon coupling in graphite 

Electron-phonon interaction theory is an extensively 
studied subject in condensed matter physics. The im- 
portance of its role is high-lighted in the theory for con- 
ventional superconductors, i.e., the BCS theory. Even 
though a general theory should be applicable to any sys- 
tem, specific and more applicable models have been de- 
veloped for metallic systems. However, a key assump- 
tion used for metallic systems, constant DOS near the 
Fermi level, is not valid for semi-metals and insulators. 
To the best of our knowledge, EPCs in semi- metals and 
insulators have not been thoroughly studied theoreti- 
cally (probably due to lack of interest). With recent de- 
velopments in graphene/graphite related researc h-^^ ^, 
EPC in semi- metals has become more important. There- 
fore, we need a model to evaluate the EPC constant in 
graphite. 

To understand the EPC in graphite, one should con- 
sider it's characteristic band structure near the Fermi 
level. Fig. 1(a) shows the crystal structure of graphite. 
Graphite has a layered structure and the stacking order 
is ABAB... In each layer, carbon atoms form strong <j 
bonds produced by sp^ hybridization while the out-of- 
plane Pz orbitals form tt bonds. Fig. 1(b) depicts the 
first Brillouin zone (BZ) and high symmetry points of 
graphite in reciprocal space. The calculated electronic 
band dispersion of graphite along the high symmetry 




FIG. 1: (a) The structure of graphite. Spheres are carbon 
atoms. Graphite shows a layered structure, which has an 
stacking order of ABAB... (b) First Brillouin zone of graphite. 
The symbols represent high symmetry points, (c) Calculated 
electronic band structure along A-L-H-A. (d) Approximated 
band structure. Each corner of hexagon is H point and z 
direction is energy. 



line, A-L-H-A, is plotted in Fig. 1(c). The band dis- 
persion within ±1 eV near the H point is almost lin- 
ear. A three dimensional view of the band dispersion is 
shown in Fig. 1(d). The point at which the two cones 
meet each other is at the Fermi energy and is called the 
Dirac point. This band structure yields a DOS which 
increases linearly with binding energy (linear DOS). Nu- 
merous studies of the band structure of graphite can 
be found both iheomiicci^^-^^-^ and experimental 
-idc3^^^^^^^^^^^^^^ 

The electron-phonon coupling theory in graphite 
should be considered within this characteristic linear 
DOS. To understand the electron-phonon coupling in 
graphite, one needs to get the real or imaginary parts 
of the self energy. If one approaches electron-phonon 
coupling through the real part of the self energy, one has 
to confront a task of finding the bare band. The bare 
band of graphite is not linear and hence much harder to 
guess in comparison with metallic systems. Even though 
it was argued that the experimentally measured band 
structure at H is linear^^, our results show that the dis- 
persion is not linear and has some parabolic character 
near the Fermi energy. Therefore, we chose to use the 
imaginary part in the analysis. Note that the real part 
of the self energy can be obtained by Hilbert transform- 
ing the imaginary part. Fig. 2(a) shows the Feynman 
diagram^^ for the lowest order EPC under considera- 
tion. Our model considers only this lowest order EPC 
in this section. One can describe the many-body effects 
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FIG. 2: (a) The lowest order Fyenman diagram for EPC. g 
is electron phonon coupling constant, k and k' are crystal 
momenta of holes, q is momentum of phonon. (b) Schematic 
diagram of the EPC process as shown in panel (a). Photo-hole 
k makes a transition to k' emitting phonon of q. Tiujq is the 
emitted phonon energy, (c) Schematic diagram for scattering 
in kz direction, (d) The imaginary part of the self energy vs. 
binding energy predicted by our qualitative theory (see the 
text). 



on quasi-particles using a self energy scheme. The imag- 
inary part of self energy is proportional to the scattering 
rate of quasi-particle. Here, we present the EPC process 
for zero (T = 0) and finite temperature cases (T ^ 0) 
separately. 



Zero temperature case 



The Hamiltonian of EPC interaction can be written 
as, 



where ^ (cfe,cr) creates (annihilates) an electron with 
spin a and momentum k while 6j {bq) creates (annihi- 
lates) a phonon u with momentum q. The scattering 
amplitude g is taken to be energy and momentum in- 
dependent, term is for phonon emission process and 
b-q^i, term is for phonon absorption. Then the imaginary 



part of the self-energy is defined as a convolution over the 
density of states: 



V 



(2) 



where V is the electronic DOS and uOj, is the energy of 
a phonon v. f and b are Fermi-Dirac and Bose-Einstein 
distributions, respectively. Since the electron-phonon in- 
teraction does not alter the spin of a conduction electron, 
spin index a is suppressed to consider only one spin di- 
rection. 

If we assume an Einstein phonon with an energy of coq 
and momentum independent coupling amplitude g, Eqn. 
(2) becomes 



^V{u + u;o)[f{uJo + cj) + 6(cJo)]) 



(3) 



For the zero temperature case, Fermi-Dirac function 
can be replaced by step function and the Bose factor is 
zero in Eqn. (3). Then, Eqn. (3) can be written as. 



S;;(c^, T = 0) = -g\[D{u - uo)e{u - uo) 

+P(cj + cjo)©(-^-^o)] 



(4) 



where is a step function, Q{x) = 0{x < 0) and Q{x) = 
l{x > 0). 

We assume a conical band structure with the Fermi 
energy at the apex of the cone. There is another con- 
ical band above the Fermi energy which is unoccupied, 
and these two conical bands form a Dirac-cone-like band 
structure as shown in Fig. 2(b). If a photo- hole with mo- 
mentum k (filled circle) is created by a photon as shown 
in Fig. 2(b), it can be filled by an electron with energy 
of uJk' = uJk — ujQ and momentum = k — q (empty 
circle) where q is the phonon momentum. The scatter- 
ing rate is proportional to the number of such k' states, 
thus the DOS at cok — ^o- Note that if the binding en- 
ergy of the photo-hole is smaller than the phonon energy 
cjo, the scattering cannot occur because there are no elec- 
trons with sufficient energy to emit a phonon with energy 
coq. Therefore, the imaginary part of self energy of photo- 
hole as a function of the binding energy is proportional to 
V{uJk—^o) and looks like a schematic shown in Fig. 2(d). 
Note that it is also possible that a phonon may scatter a 
photo- hole in the kz direction as shown in Fig. 2(c). The 
outcome is not much affected by the c-axis scattering due 
to the weak dispersion of the tt band along kz direction. 

Once is obtained, one can obtain the real part of 
self energy by Hilbert transforming T^" . The electron- 
phonon coupling parameter, A is defined as. 



A: 



duo 



(5) 



a;=0 



At the K point, the bonding and anti-bonding bands are 
split because of the inter-layer interaction of graphite^!. 



Considering the small inter-band scattering of the photo- 
hole by a phonon between bonding and non-bonding 
bands at K point, the above self energy can be extended 
to the double band case at K. This double band case 
was investigated in our previous work^^. 

We also note that is not affected seriously by the 
detailed shape of near cc; = because increases 
linearly. This aspect was considered in calculating 
for K and H points^. It was argued that H^p is some- 
what different at K and H because the band structure 
at K is parabolic near the Fermi level while that at H is 
linear. Meanwhile some difference between K and H cer- 
tainly exists that affects the detailed shape of E^p near 
the Fermi energy, the effect on the EPC constant A should 
to be negligible because the contribution comes mostly 
from the high binding energy side. 



2. Finite temperature case 

We now move onto the finite temperature case. In the 
case of metals, there is an easy way to extract the EPC 
constant A from temperature dependent data through a 
simple formula^^^^. The formula is derived under the as- 
sumption that the electronic DOS near the Fermi energy 
is constant, which is not the case for graphite. Here, we 
investigate the temperature dependence of S^p theoret- 
ically to determine if one can easily extract A from the 
temperature dependence data. It turns out that a simple 
formula such as the one for metals can not be formulated. 
However, we show some possibility of estimating EPCs 
from the temperature dependent data. 

The imaginary part of the self-energy by electron- 
phonon coupling at finite temperature was shown in Eqn. 
(3) of the previous section. Note that for the high phonon 
frequency (for example, A'-^ or E2g mode in graphite) 
the Bose factors can be neglected for the temperature 
range over which we performed our experiments (lOK- 
225K). We consider not only these high energy phonons 
but also low energy phonons. Unfortunately, the temper- 
ature dependence of S^p in Eqn. (7) cannot be reduced 
to a simple form as the one for a metal^ and extracting 
the A from the temperature dependence of S^p is not an 
straight forward task. However, one can still obtain infor- 
mation from the temperature dependent data. If the low 
frequency phonon mode participates in EPC, the tem- 
perature dependence of S^p near the Fermi level should 
be strong while the high frequency phonon modes should 
contribute little to the temperature dependence. In ad- 
dtion, the Bose factor in Eqn. (3) is not negligible and 
Sgp will show clear difference at different temperature. 
Therefore, once the g value is known, one can roughly 
identify which phonon mode contributes the most to EPC 
by fitting the temperature dependent data. 
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FIG. 3: (a) The Feynman diagram for the lowest order for 
electron-hole pair creation process. Photo-hole k makes a 
transition to a hole k' ^ and an electron ke and a hole ku 
are created conserving momenta and energy, (b) Schematic 
diagram for electron-hole pair creation in phase space, (c) 
Schematic diagram for electron-hole pair creation considering 
the scattering along kz direction, (d) Hatched area repre- 
sents the possible A/c and Aa; for electron-hole pair creation. 
Dashed line represents the possible A/c and Acj for photo-hole 
transition. These two area slightly touch each other at line, 
therefore there is no available phase space. 



B. Electron-electron interaction in graphite 

Electron-electron interactions can also affect the quasi- 
particle lifetime or the imaginary part of self energy. We 
consider the lowest order scattering in electron-electron 
interaction occurring via electron-hole pair creation. The 
Feynman diagram for this scattering process is shown in 
Fig. 3(a). When the photo- hole, of crystal momentum is 
k, is created, the hole makes the transition to k' by cre- 
ating another hole, kh, and electron, kg. Fig. 3(b) shows 
this electron-hole pair creation process in the E — k phase 
space. From the figure, one can see that the electron-hole 
pair creation process is negligible under the linear DOS 
of graphite near the Fermi energy. Let us consider the 
available phase space for electron-hole pair creation in 
Fig. 3(b). If we plot the energy difference (Acj = uok' —^k 
where cok' and cok are the energies of the holes with k' and 
k, respectively) as a function of momentum difference 
(Ak = k' — k), the possible transitions occupy the area 
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below the dashed hne in Fig 3(d). In a similar way, one 
can find that electron-hole creation process occupies the 
hatched area in Fig. 3(d). Photo- hole decay through the 
electron-hole pair creation can occur only when the two 
conditions are met, that is, where the phase spaces for the 
two processes overlap. They overlap only on the dashed 
line as can be seen in Fig. 3(d). Therefore, the available 
phase space volume for the decay through electron-hole 
pair creation is zero. Note that, this is true only near the 
Fermi energy where the band structure can be approx- 
imated by Dirac cones. For the photo- holes at higher 
binding energies, the available phase space volume is be- 
comes non-zero. This fact was previously pointed out by 
Moos et al.^ If we limit our discussion to the low energy 
dynamics in graphite, the effect of the electron-electron 
interaction can be neglected. 



C. Other scattering mechanisms 

There are other mechanisms in graphite that may 
contribute to the quasi-particle scattering such as plas- 
mons, impurities and defects. Xu et al. suggested 
that plasmons may be the main source for the quasi- 
particle scattering in graphite^. However, Spataru et 
al. showed that elctron-hole pair creation should be a 
more dominant mechanism than electron-plasmon inter- 
actions for electron scattering in graphite^. Since we 
have shown in our earlier work^"^ that electron-phonon 
interaction is more dominant than electron-hole pair cre- 
ation based on a phase space argument, we may conclude 
that electron-plasmon interaction is much weaker than 
electron-phonon coupling and thus may be neglected. 
Impurity and defect scattering can also contribute to the 
scattering rate in graphite. These scattering mechanisms 
also have a rate that is proportional to the electronic 
DOS as in the electron-phonon coupling case and thus 
increase the slope for the imaginary part of the self en- 
ergy. This fact tells us that, if one wants to study the 
electron-phonon coupling, the experiment should be con- 
ducted on clean single crystalline graphite. In our case, 
we used natural graphite single crystals which have supe- 
rior quality to crystals used in other experiments. As a 
results, we did not observe any defect related states^^^^^ 
and we therefore believe that defect or impurity scatter- 
ing is minimal. 

In short electron-phonon coupling should be mecha- 
nism, other mechanisms are suppressed due to lack of 
phase space (electron-hole pair) or high quality of the 
crystal (low impurity /defect levels). 

III. EXPERIMENT 

ARPES experiments were performed at Beamline 7.0.1 
of the Advanced Light Source. We used very high quality 
natural graphite single crystals with sizes larger than ^ 
1 cm. Samples were cleaved repeatedly ex situ by tap- 



ing method until a fiat surface without large flakes were 
obtained. Samples were subsequently introduced to the 
ultra high vacuum chamber and annealed at 900° C for 
30 minutes in a vacuum better than 6.0x10"^^ Torr. to 
clean the surface. The energy resolution was ~ 40 meV. 
The chamber pressure was better than 5.0x10"^^ Torr. 
during the measurements. We found that typical size of 
the flat regions without flakes was smaller than 200 /am. 
Therefore, we exploited the small beam spot 50 fim) 
to probe flat region. 

We took kz = H data at 20K with a photon energy 
of 103.4 eV to obtain the electron-phonon coupling by 
analyzing the peak width as a function of the binding 
energy. This is essentially the same as what we reported 
earlier^^ but at the H point. In addition, we performed 
temperature dependent experiment at the K point with a 
photon energy of 85 eV. ARPES data was taken at 25K, 
75K, 125K, 175K, 225K. We started measuring at 225K 
and lowered the temperature. After having measured at 
25K, we annealed the sample again for ^ 30 seconds, at ~ 
900° C and measured again. Comparison of the data be- 
fore and after annealing showed essentially no difference, 
indicating there was no surface contamination during the 
measurement. For comparison, graphene data were also 
taken at the K point. The graphene sample was epitax- 
ially grown on 6II-SiC in situ as reported elsewhere^S. 
Electronic band structure calculation was done by using 
the SIESTA code based on pseudo-potential method. 



IV. RESULTS AND DISCUSSION 

A. Low temperature case 

Fig. 4(a) shows measured ARPES spectral function 
along the L-H-A symmetry line. The non-bonding band 
(NB) and bonding band (BB) are degenerate at the H 
point whereas they are split at K point. We could iden- 
tify only one peak in the energy (EDCs) and momentum 
(MDCs) distribution curves. We also took data with dif- 
ferent photon energies, to insure that we were really at 
the H point. The electronic band near the Fermi energy 
shows a linear dispersion as predicted in the band calcu- 
lation in Fig l.(c). However, we also note that the band 
dispersion very near the Fermi level shows some parabolic 
component contrary to what is expected from the theory. 
This could be due to kz broadening caused by the finite 
escape depth. We also note that there are no evidence 
defect-induced states that were reported earlier^^. This 
indicates that our natural graphite single crystals are of 
very high quality. Almost, the negligible background of 
our data even at high binding energy further supports 
the high quality of our sample. This means that defect or 
impurity contribution to the scattering rate is very small 
and we may only consider the electron-phonon coupling 
effect. 

Fig. 4(b) shows the EDC from the k point indicated 
by the arrow in panel (a). The line shape of the EDC 
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FIG. 4: (a) ARPES data taken at H point of graphite along L- 
H-A direction, (b) The EDO at the k point marked the arrow 
in panel (a). Circles are the experimental data and thick line 
is the fit for which finite escape depth effect in photoemission 
process has been considered (see the text), (c) Contributions 
from different kz points due to the finite escape depth effect, 
(d) The model fitting function with finite escape effect consid- 
ered but without lifetime effect, (e) ARPES data at K from 
epitaxially grown graphene on 6H-SiC. (f) The EDC (at the 
k point marked by the arrow in panel (e)) shows symmetric 
lineshape unlike that from graphite. The EDCs can be fit- 
ted with a single Lorentzian and constant background (thick 
line) . 



is very asymmetric. As was the case for the K data^^, 
we can understand this asymmetry as follows. Though 
we tuned the photon energy to probe the H point of 
graphite, the finite escape depth of the photoelectron 
yields an uncertainty in kz , A/c^ = 1/m where /i is the es- 
cape depth. Therefore, there is contribution from other 
kz values which is illustrated in Fig. 4(c). As the BB 
and NB have finite kz dispersions, the contribution from 
other kz values results in broadening of the spectral func- 
tion. The fact that the BB has more kz dispersion gives 
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FIG. 5: (a) Extracted HWHM as a function of the binding 
energy for NB (triangles) and BB (filled triangles), (b) Calcu- 
lated pDOS for NB and BB (c) Constructed T." from the data 
in panel (a) for the low energy region and pDOS in panel (b). 
pDOS is scaled so that it matches the experimental HWHM 
at 0.9 eV. (d) using Hilbert transform of T." . The EPC 
parameter is ~ 0.23. 



more broadening on the higher binding energy side as is 
seen in Fig. 4(b). 

Fig. 4(d) depicts a model spectral function when all 
these effects are accounted for. Only when such effects 
are considered, can one extract the true lifetime broad- 
ening. We used /i = 7A for the fitting^, and the model 
function in Fig. 4(d) is convolved with a Fermi function 
and a Voigt function with the Gaussian width set to the 
total energy resolution of 40 meV. In addition, we in- 
troduced binding energy dependent Lorentzian width for 
the Voigt function considering the observation from the 
K data that the Lorentzian width linearly increases as a 
function of binding energy 

Even though this kz uncertainty is a general prop- 
erty of ARPES measurement, perfect 2D material such 
as graphene should not show this escape depth effect in 
their ARPES data because it has no dispersion in the kz 
direction. To ensure that the asymmetric line shape in 
graphite is indeed from the finite escape depth effect, we 
took graphene ARPES data and check if the line shape 
is symmetric as expected. Fig. 4(e) shows ARPES data 
along the M-K-T direction of graphene which was epi- 
taxially grown on 6H-SiC. Fig. 4(f) is an EDC curve 
from the k point arrow marked in panel (e). The EDC 
shows very symmetric line shape contrary to the EDC in 
panel (b). One can fit this curve with a single Lorentzian 
with constant background as shown with the thick line 
in Fig. 4(f). An almost perfect fit strongly supports the 
idea that the asymmetric line shape of graphite data is 
indeed from the finite escape depth effect. 

Fig. 5(a) shows the half width at half maximum 



7 



(HWHM) found by fitting our model to tiie EDCs along 
the high symmetry line H-A. Filled and empty sym- 
bols represent BB and NB widths, respectively. There 
is almost no difference between the BB and NB widths. 
The width increases linearly as a function of the bind- 
ing energy. We find that the width shows no high order 
dependence such as ~ o;^. This also indicates that the 
EPC is the dominant decay channel in graphite as ex- 
pected from our model. Yet, observation of very weak or 
no kinky feature at the optical phonon energy of 0.2 eV 
shows that EPC is very weak in graphite. On the other 
hand, the width converges to zero as the binding energy 
goes to zero, which means that momentum mixing due 
to impurity or defects is minimal, supporting again the 
high quality of the samples. 

One can extract the EPC parameter from the deriva- 
tive of Ti' at 00=0. Conventionally, one obtains the T^' 
from the difference between the experimental dispersion 
and the bare band. In graphite, this is a difficult task 
because the bare band may not be linear. On the other 
hand, even though harder, one can get by Hilbert 
transforming . To do the Hilbert transformation, 
we need to know T^" over the entire energy range. As 
this is not the case, we use scaled partial electronic DOS 
(pDOS) as E^', assuming that is approximately pro- 
portional to pDOS^"^. Fig. 5(b) shows pDOS of NB and 
BB. Fig. 5(c) is the T^" for NB band, obtained from the 
experimental data and calculated pDOS for NB. Hilbert 
transform of it gives the shown in Fig. 5(d). Accord- 
ing to Eqn. (5), we can find the electron phonon coupling 
parameter from the energy derivative of at a; = 0. The 
resulting value is A ^ 0.23, very similar to the value of 
A = 0.2 for the K-T direction reported in our previous 
worki^. This value is larger than the calculated value of 
0.075 for graphene^ but much smaller than the previously 
reported value for graphite^^. In addition, this value is 
consistent with the value of 0.21 calculated with a rea- 
sonable scattering amplitude g^. Therefore, we conclude 
that EPC constant A is also small at the H point. 



B. Finite temperature case 

Fig. 6 shows temperature dependence ARPES data at 
K point of graphite, which were taken at 225, 175, 125, 
75, and 25K. One can clearly distinguish the NB from the 
BB in each panel. Note that the binding energy differ- 
ence between the NB and the BB is about 0.8 eV. Overall, 
the data do not appear to show much temperature de- 
pendence. To see this quantitatively, we performed the 
same lineshape analysis we developed on the data. Every 
EDC from -0.2 to A of each panel in Fig. 6 is fitted 
with our model function and HWHM is extracted. 

Extracted HWHM vs. binding energy at different tem- 
perature is plotted in Fig. 7. Overall, HWHMs linearly 
increase proportional to binding energy. All HWHMs 
are quite similar to each other and one can safely say 
that there is no clear temperature evolution of spectral 



function. This already indicates that the energy of the 
phonon mode which is involved in electron-phonon cou- 
pling in graphite is very high compared to the tempera- 
ture scale of our measurement 225K. 

Calculated imaginary part of self-energy is used to fit 
the extracted HWHM of the 25K data (Fig. 8(a)). In 
fitting the data, we assumed an Einstein phonon of coq 
= 200 meV. The imaginary part is supplemented with a 
constant plus an energy-dependent term See = in 
order to simulate the electron-electron interactions. The 
partial DOS of the NB band were calculated from the 
LDA calculation. The coupling amplitude ^ is a fitting 
parameter. We find that g is ?^0.39 eV and electron- 
electron interaction pre-factor is ^0.004. Note that neg- 
ligible electron-electron interaction near the Fermi energy 
is confirmed as predicted earlier in our model. 

With the fit result, one may try to evaluate the 
electron-phonon coupling constant A from the g value. 
By Hilbert transforming the model fit function of S^', we 
obtain T^' . Derivative of at a; = as in Eqn. (4) gives 
us A of ^0.14. This value is smaller than the value of 
A = 0.2 at i^^ which was obtained from the experimental 
data^^. A key difference between the two methods is that 
while we assume an Einstein phonon of a; = 200 meV, 
no such assumption is used in transforming the exper- 
imental data. However, the experimental data is more 
susceptible to systematic errors, especially at very low 
binding energy range (where the line shape is affected by 
Fermi function). Since the low energy range has more 
effect on A and a theoretical result shows u = 200 meV 
is the dominant phonon^, A = 0.2 probably gives us the 
upper bound for the electron-phonon coupling. 

Other panels in Fig. 8 show HWHMs and fitted model 
function at different temperatures. The fitting was con- 
ducted on NB bands along M-K near the Fermi energy 
at each temperature. As could already been seen in Fig. 
7, fitting of the HWHM results in negligible temperature 
dependence. This indicates that the phonon involved in 
the coupling has much higher energy scale than 225K. In 
fact, the 200 meV bond-stretching mode may be the most 
dominant one as the ab initio calculation on graphene 
shows 



V. CONCLUSION 

We present high resolution ARPES data taken at the 
H point of natural graphite single crystals. The graphite 
bands shows a linear dispersion as predicted in LDA cal- 
culation and the NB and BB are degenerate. First, we 
considered various scattering mechanisms in graphite. 
We deduced theoretical formula for the scattering rate 
by phonons. We find that the scattering rate by EPC in- 
creases linearly with binding energy due to the linear den- 
sity of states. Electron-electron interactions in graphite 
are negligible in the low binding energy region where the 
band dispersion is linear. The impurity or defect scatter- 
ing rate is also proportional to binding energy because of 
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FIG. 6: ARPES data taken along M — K — F Sit different temperatures. (a),(b),(c),(d), and (e) were taken at 225K, 175K, 
125K, 75K, and 25K, respectively. 
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linear DOS. We show that all effects other than phonon 
scattering are negligible. Second, with the finite escape 
depth effect in photoemission process considered, we ex- 
tracted T^" from the EDCs of the NB and BB separately. 
Finally, we approximated Tj" by combining the experi- 
mental HWHM Y." and calculated partial DOS. The ob- 
tained T^" is converted to Y^' through a Hilbert transform. 
The extracted EPC parameter at is ~ 0.23, which is 
small, consistent with the value 0.2 from the K point in 
our previous work. This small EPC parameter is also 
consistent with very weak kinky features in our data. 

In addition, we conducted temperature dependent 
ARPES measurements on the graphite K point. The 
temperature dependent data shows no notable evolution 
in the EDO lineshape within the temperature range (25- 
225K). Analyzing the experimental temperature depen- 
dence of peak widths and simulated temperature depen- 
dence, we conclude that the dominant phonon mode in 
EPC in graphite is much larger than the temperature 
scale of our experiment 225K. This is consistent with the 
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FIG. 8: HWHMs extracted from ARPES data at different 
temperatures are fitted by calculated imaginary part of self- 
energy. Panel (a),(b),(c) and (d) are for temperatures of 25, 
75, 125 and 175K, respectively. The solid line in each panel 
is the best fit to the experimental data by calculated imag- 
inary part of self-energy. The calculated imaginary part of 
self-energy includes electron-phonon and electron-electron in- 
teraction terms. The electron-phonon coupling constant and 
electron-electron interaction pre-factor were used as fitting 
parameters (see the text). 



notion that the phonon mode in electron-phonon cou- 
pling in graphite is the 200 meV optical phonon mode as 
is the case for graphene. 

Even though electron-phonon coupling has been heav- 
ily studied, most of these studies were focused on metallic 
systems where the density of states near the fermi level 
is approximately constant. Such is not generally true, es- 
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pecially for semi- metals. The formulas discussed in this 
work are very general and can be used for any shape of 
electronic density of states. It should therefore be useful 
in the future studies on semi- metals. 
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